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Abstract: We study the background cosmological dynamics with a three component 
source content: a radiation fluid, a barotropic fluid to mimic the matter sector and a 
single scalar held which can act as dark energy giving rise to the late-time accelerated 
phase. Using the well-known dimensionless variables, we cast the dynamical equations into 
an autonomous system of ordinary differential equations (ASODE), which are studied by 
computing the hxed points and the conditions for their stability. The matter huid and the 
scalar held are taken to be uncoupled at hrst and later, we consider a coupling between the 
two of the form Q = y^2/3Kl3pm^ where pm is the barotropic huid density. The key point of 
our analysis is that for the closure of ASODE, we only demand that the jerk, T = VV"/V'"^ 
is a function of acceleration, 2 ; = —MpV'/V, that is, T = 1 -|- f{z). In this way, we are 
able to accommodate a large class of potentials that goes beyond the simple exponential 
potentials. The analysis is completely generic and independent of the form of the potential 
for the scalar held. As an illustration and conhrmation of the analysis, we consider f{z) 
of the forms {p ~ {p ~ to numerically compute the evolution of 

cosmological parameters with and without coupling. Implications of the approach and the 
results are discussed. 
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1 Introduction 

The evidence of dark energy being present in the cosmological soup comes from various 
observations such as the distance measurements from type la supernovae [1], baryon acoustic 
oscillations [2], measures of the Hubble constant [3] and the precision measurements of 
cosmic microwave background [4]. The data from these observations consistently favours 
the spatially flat ACDM model with a currently accelerating phase. Here, the cosmological 
constant A is used as an ingredient responsible for dark energy with the equation of state 
parameter rcA = — 1 to produce the late-time accelerating phase. However, an exceptionally 
small value of the cosmological constant with AL^ ~ (fine-tuning) and the relatively 

equal order of magnitude values of Ha ~ 0.7 and ~ 0.3 (coincidence) suggest that dark 
energy could be more than just having a cosmological constant in general relativity^. 

There are many possibilities that have been explored to account for the late-time ac¬ 
celeration of the universe (for a non-exhaustive sample see refs. [6-17]). A popular and 
simplest approach is using a light scalar field termed “quintessence” as a candidate for dark 
energy. In the quintessence models, the scalar field is either taken to be slowly rolling along 
a shallow potential (“thawing” class) [18] or else with a “scaling” behaviour [19-22] (which 
contains “tracker” as a special case) in which it tracks the dominant fluid component and in 
the later stage exits the scaling regime with a quasi-constant equation of state parameter. 
The form of the self-interacting potentials for the scalar field is usually phenomenological 

^It should be noted that apart from these concerns, there are no practical issues concerning the grav¬ 
itational effects and cosmological implications of having a cosmological constant. See ref. [5] for a field 
theoretic viewpoint of the cosmological constant problem. 
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or else is motivated as arising from the theories with extra degrees of freedom such as 
Kaluza-Klein reduction of a higher dimensional theory [23, 24] or from string theoretic [25] 
considerations. Further we can also consider, in principle, an interaction between dark en¬ 
ergy (the scalar field here) and matter fluid (including dark matter). Such an interaction is 
again either introduced phenomenologically to address the coincidence problem (also affect¬ 
ing structure formation in new ways) or can be motivated by considering modified gravity^ 
models in the Einstein frame [26-31]. The couplings that arise from such considerations 
affect the matter sector with an effective metric of the form [32-35]: 

9ah = C{(l))gah + ( 1 - 1 ) 

where and D{(j)) refer to the conformal and the disformal part of the coupling respec¬ 
tively. With these model inputs, we can study the background and perturbed cosmological 
dynamics for parameter estimation as well as to infer the concordance with observations. 
The background evolution is usually and effectively studied by transforming the cosmologi¬ 
cal dynamical equations to an equivalent non-linear autonomous set of differential equations 
[36]. The only tractable analytical approach is to look at the fixed points of the equivalent 
system and their stability conditions. These fixed points can be associated with different 
cosmological phases which can be either early-time (source), transients or late-time (sink) 
solutions. The numerical integration of the differential equations with given initial condi¬ 
tions reveals the phase space dynamics and the evolution of cosmological parameters. We 
can compute the background quantities such as the evolution of Hubble parameter, equa¬ 
tion of state parameter, luminosity distances etc. and compare these with observations. 
The perturbed evolution requires running modified versions of popular CMB Boltzmann 
codes. 

However, many subtleties hide in the modelling of the dark sector and one requires a 
careful analysis of each step one at a time which brings us to the motivation (for yet another 
paper) and the specific point of view adopted here. In this work, 

• We study the background cosmological evolution in a fairly generic setting. The 
analysis encompasses a large class of potentials and the interacting dark sector models 
including those arising out of the modification of gravity sector all (except disformal) 
in one go. 

• We consider a three component system with radiation, a barotropic fluid and the 
quintessence scalar field with a (quite) generic potential. This is as close as we can 
get to the actual composition of our universe that can be handled analytically. 

• We present a classification of a large number of potentials for the quintessence field 
beyond the simple exponential potential. There exist three main classes for all such 
potentials based on the properties of the function T — 1 = VV" — 1. 

^In some circles, the presence of dark energy is synonymous with modifying the gravity sector rather 
than having an extra source in Einstein gravity. However, conformal transformations can bring most of the 
modified gravity actions back to the canonical Einstein form with extra source degrees of freedom. It is then 
a matter of debate to ask which one of these conformally related frames is more physical or fundamental. 
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We then proceed to understand the cosmological dynamics through dynamical system 
analysis in terms of computing the fixed points, understanding their properties and 
stability conditions. 


• Finally, we resort to numerical techniques and show the evolution of cosmological 
parameters for a few examples in which tracking or thawing behaviour is observed. 

Most of the previous studies on similar lines [37, 38] consider only two components - the 
scalar field and a matter fluid - focusing only on the late-time accelerated regime where the 
contribution of radiation is negligible. However even though it is subdominant at present, 
in order to obtain a correct tracking regime [39] from early epochs, we have to include 
in radiation from the start. We have not seen a complete and generic analysis of three 
component system anywhere in the literature and it is therefore interesting to see what we 
can infer from it. 

The outline of the paper is as follows: We introduce the machinery of an equivalent 
dynamical system in Section 2.1, classification of the potentials in Section 2.2 and briefly 
review the linear stability analysis in Section 2.3. In Section 3, we compute the fixed 
points, the conditions of their existence and stability and cosmological implications with 
and without coupling. Using numerical techniques in Section 4, we consider some examples 
to demonstrate and confirm the analytical results obtained in the previous section before 
summarizing the conclusions in Section 5. 

2 Theoretical Setup 

In this section, we introduce the theoretical framework of the (interacting) quintessence 
models with a brief review of the dynamical systems analysis relevant for our purposes 
(see ref. [36] for an excellent introduction). We, first, develop the machinery to go from 
the cosmological dynamical equations to an equivalent system of dimensionless phase space 
variables. These satisfy an autonomous set of differential equations (ASODE). We then 
mention the rules involved in the stability analysis of a non-linear autonomous system. 

2.1 Prom cosmological dynamical equations to ASODE 

We assume a spatially flat, homogenous and isotropic background which gives the diagonal 
Friedmann metric, gab = ( —1, a^, a^, a^). Assuming the framework of general relativity and 
the sources, 

• a scalar field 0 , minimally coupled to gravity: = </>^/2 — V{(f>) with the equation 

of state, = (70 - l)p^ where p^ = 0^/2 + V{4>) and = <^2/2 - !/(</)), 

• a barotropic fluid (matter sector) with pm = WmPm = (7 ~ ^)pm and 1 < 7 < 2 , 

• a radiation fluid for which Pj. = Pr/3, 

the dynamical equation for this metric is given by. 



( 2 . 1 ) 
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where we have the standard dehnitions, H = a/a and Mp = l/(8'7rG). With this, the 
background evolution is also required to satisfy, 

as the flatness constraint. Further, the stress-energy tensor for the sources has to be 
divergence-less. In the absence of any coupling between the components, the stress-energy 
tensors for the components are divergence free individually. However, in the presence of 
any coupling, the divergence-free condition holds in an added manner. This leads to the 
following set of (continuity) equations: 


P0 -h 3H'y^p^ = Q 

(2.3) 

Pm T Pjyi — Q 

(2.4) 

Pr + 4Hpr = 0, 

(2.5) 


where Q denotes the coupling between the scalar held and the matter sector. One can 
construct the coupling purely on the basis that at least Q should depend on the energy 
densities and other covariant quantities. As remarked before, we can also get a coupling 
between the scalar held and the matter sector by writing a modihed (scalar-tensor) gravity 
theory in canonical Einstein form (see refs. [26, 27]). Locally, the modihed continuity 
equations imply the rate of energy exchange between the two components which can be 
classihed as 


Q 


> 0 
< 0 


dark matter —)• dark energy 
dark energy —>■ dark matter 


( 2 . 6 ) 


The dynamical equation for the scalar held with this coupling is given by. 


+ + = Q//) 


(2.7) 


obtained by substituting for in Eq. (2.3) and hence is not independent of the continuity 
equation. This completes the set of dynamical equations for our system. For our purpose 
of the tracking the behaviour of and r/;^, it is convenient to introduce the dimensionless, 
phase space variables [40, 41]: 

2 — . 2 — ^ . 2 _ o — Pr _ fr, Q\ 

^ ~ ^ “ 3H^Ml ^ ^ 

such that the dimensionless energy density of the held (j) is given by = x'^ + and 
Wfj) = (x^ — 2 /^)/(x^ + y^)- The Friedmann constraint can be written as 

= 1 — (2.9) 

This constrains the physical phase space with the conditions, 

0 < < I ^ 0 < x'^ + y'^ + < 1. 


- 4 - 


( 2 . 10 ) 







Besides, since 0 < < 1, we have 0 < < 1 as well as |x| < 1. Additionally, with 

positive cosmic expansion, that is, // > 0, we have y > 0. Also, 0 < < 1 gives |tt| < 1 

but we shall see that the system of equations will be symmetric under u —)• — u and hence 
the physical phase space is defined by, 

:= {(x, y, ti) ; 0 < < 1,0 < < 1, |x| < 1, y > 0, ti > O} (2.11) 


which is a bounded upper semi-disk. This is enough for the potentials that are constant or 
belong to simple exponential class, For the potentials beyond these simple 

classes, we need to introduce more variables. The first such variable [42, 43] is, 

z = -Mpy (2.12) 

which vanishes for a constant potential and is constant for an exponential potential. In 
terms of these set of variables {x, y, u, z}, we can recast the dynamical equations as: 


dx 

m 

dy 

dN 

du 

dN 

dz 

m 


-3x + v -y z + -X 


,21 


(2 - 7 )x -k 7(1 - y ) (4 - 37 ) — 


+ 


Qx 


/3 3 

-^-^xyz+-y 


,21 


(2 - 7 )x -h 7(1 - y ) (4 - 37 ) — 


-2u + -u 


, u 


(2 - 7 )x -h 7(1 - y ) (4 - 37 )-^ 


—V^xz^iv — 1)] r = 


vv" 

yi2 


(2.13) 

(2.14) 

(2.15) 

(2.16) 


where we have now introduced N = J da/a as the time parameter which is also a mono- 
tonically increasing function. For the closure of the above set of equations, we now demand 
that F — 1 = f{z) where f{z) is some arbitrary function of 2 . Otherwise, one needs to 
consider an additional ODE: F^tv = ... and the phase space becomes more complicated. 
However, even within this constraint of the jerk, F being some function of the acceleration 
z, we are able to accommodate a very large class of potentials beyond a simple exponential 
class. With this framework, we now wish to understand the dynamics of the variables, 
{x,y,n, z} as per the above autonomous set of differential equations^(ASODE). 


2.2 Classification of potentials 

With the identification of (F — 1) with f{z) in Eq. (2.16), we still have a large class of 
potentials at hand for the scalar field. The mapping of the potentials V{(/) with the function 
f{z) and the form of the right hand side of Eq. (2.16) i.e., dz/dN = —y/6xz‘^f{z) provides 
a way to classify them as follows: 


• Potentials for which /(z* = 0) = 0, e.g., f{z) = yz^ or (sinh(/r 2 ;))" with a > 0, 

• Potential for which /(z^,) = 0 with z* ^ 0. In this case, we can have two further 
classifications. 

®The system, x'i = fi{xj) forms an autonomous set of differential equations (ASODE) when the functions 
fi{xj) do not depend on the evolution parameter explicitly. 
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— z‘^f{z)\z=o / 0 e.g., f{z) = (/i — z^)lz^. Then x = 0 is necessary for an allowed 
fixed point. 

- z^f{z)\z=o = 0 e.g., f{z) = in- or (/x - z°^)/z. 

• Potentials for which z* is not finite e.g., f{z) = iijz'^ with 0 < a < 2. 


We can get the potentials connected with an f(z) easily by computing V{z) and then 
substituting for z{4>). From the definition z = —MpV'/V it is straightforward to see that. 



a. 

II 

(2.17) 

which gives 


(2.18) 

and similarly we have 


(2.19) 

which solves to 

(M 

1 

II 

(2.20) 


In principle, we can now invert the above expression to get z{(f)) which can be substituted 
in Eq. (2.18) giving V{(j)). However, it is important to note that the mapping between the 
potential and f{z) is not completely unique. For example, in the case of V{(p) = Vbe“® ^ 
[43], we have f{z) = j^jz which contains no information about the parameter a in the 
potential. Thus, a does not affect the dynamics even though the potential has a very 
different behaviour for a < 0 and a > 0 with a given j3. The dynamics is only sensitive to 
f{z) and any parameter estimation run using the background level observational aspects will 
not determine the potential completely. Hence, we shall just focus on f{z) in our analysis 
which in its simplest form can be written as f{z) = fj,{a + bz‘^) where the parameters 
{fj,, a, b, c} are real. 

2.3 Linear stability analysis of an autonomous system 

We can assess the dynamics of the system in the phase space by looking at the fixed or 
critical points for the ASODE and analysing them for their stability at the linear level. 
Essentially, this means that for a set of equations, x' = fi{xj), we find the set of points, 
{(x^'^^)} for which all vanish simultaneously. The next step is to consider a small 

perturbation in the neighbourhood of each fixed point, that is, taking x) ' + 5xi and ex¬ 
panding fiix^j^'^ -|- 6xj) to linear order thus linearizing our equations in order to judge the 
stability of the given equilibrium point. Then, depending on whether the linear pertur¬ 
bations decay(grow) with time N or decay in one direction while grow in the other, the 
critical point can be a future (past) attractor or a saddle point. This inference essentially 
follows from the study of the eigenvalues Ai of the Jacobian, 
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• If Re(Aj) 7 ^ 0 Vi, the point is termed as a hyperbolic critical point. 

— With Aj G Reals and 

(a) all Xi < 0, the critical point is a stable node (future attractor) or a sink. 

(b) all Xi > 0, the critical point is an unstable node (past attractor) or a source. 

(c) Aj < 0 < Xj for some i,j the critical point is a saddle point (marginally 
stable), that is, a source in one direction and sink in the other. 

— If eigenvalues are complex numbers with A, A* = u ±iuj and oj ^ 0, then if 

(a) < 0, the critical point is a stable spiral with damped oscillations. 

(b) ly > 0, the critical point is an unstable spiral with growing oscillations. 

• If Re(Ai) = 0 for some i, the point is a non-hyperbolic critical point. In this case 
if the imaginary part is finite, it is associated with free oscillations. If such a point 
has at least one positive eigenvalue, it is an unstable point. Otherwise, one has to go 
beyond the linear stability theory and use the centre manifold theory in order to infer 
the stability criterion. 

By knowing the critical points and their nature in the phase space corresponding to a given 
cosmological model, we can correlate the attractors with cosmological solutions that can 
decide the end point and/or the origin of the cosmic evolution. The meta-stable saddle 
equilibrium points which are not local extrema can be associated with equally important 
transient cosmological solutions. Thus, in a convenient trade-off, instead of studying the 
cosmological dynamics given by the Hubble parameter, H{t) and the densities of the sources 
Pi{t), we opt to study the equilibrium points of an equivalent ASODE in the phase space 
of the cosmological model. 

3 Fixed point analysis and cosmological implications 

We now turn to the computation of fixed points, their existence and the stability conditions. 
We first look at the case without any coupling between matter and the scalar field after 
which we look the effects with the coupling turned on. We also give the cosmological 
implications of each of the fixed points which can be associated with source/sink/transient 
cosmological behaviour according to their stability. 

3.1 Dynamics without any interaction between the components 

Consider the case when there is no interaction between matter and the scalar field, that 
is, taking Q = 0 in Eq. (2.13). It is clear from Eq. (2.16) that for z to be stationary, we 
need either of x, z or f{z) to vanish. Hence we can look at the fixed points case by case, 
that is, when (i) z = 0, (ii) z = z* for which /(z*) = 0 and (iii) z = Za (arbitrary) which 
requires x = 0. We find that there are 12 fixed points in this system. These along with 
their eigenvalues in the linearised system are tabulated in Table 1. The conditions for their 
existence, nature (stability) and the cosmological parameters rCeff, w^j, and at those 
points are given in Table 2. 
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Points 

X 

V 

U 

z 

Ai 

A2 

A3 

A4 

Pi 

0 

1 

0 

0 

-2 

-37 

-l±yi-lz^fiz)\z =0 

p± 

±1 

0 

0 

0 

3 

1 

Kl 

3a 

P3 

0 

0 

0 

Za 

0 

37 

2 

3 a 

2 

b 

2 

P4 

0 

0 

1 

Za 

0 

2 

-1 

b 

P5 

0 

0 

0 

Z* 

0 

37 

2 

3 a 

2 

b 

2 

Pe 

0 

0 

1 

Z* 

0 

2 

-1 

b 

P± 

±1 

0 

0 

Z* 

1 

3a 

6 t\/ 62 * 

2 

^VQz^df^ 

Ps 

Vg 

A-l 

0 

Z* 

(^2-4) 

2 

- 37 

(^*-6) 

2 

1 

* CO 

Pg 

r~8~ 
V ^ 

2 

'/iZt 

1 

> 


6 

1 _L V 64-1522 

2 2 zt 

—Az^^df^ 

Pio 

\/67 

2 Z^: 

\/6o7 

2 zt 

0 


b 

2 

3a 1 3 \/“(2472-^2(9^_2)) 
4^4 2* 

—Sjz^df^, 


Table 1: Fixed points and the corresponding eigenvalues for the i?i)-niatter-radiation system 
without coupling, i.e., Q = 0. Here 0 = 2 — 7 and 6 = 4 — 87. 

Properties of the fixed points The twelve fixed points are divided according to the 
classification of the potentials in Section 2.2. The points Pi _4 exist for all type of potentials 
while the points P 5-10 exist for only those potentials for which f(z) vanishes for some finite 
z*. The properties of these fixed points with their cosmological implications are as follows. 

Pi This point is the stable attractor solution when z'^f{z)\z=o > 0 (stable spiral for 
z'^ f{z)\z=o > 3/4) and is also associated with the accelerated de Sitter phase with rcefi = —1 
and dominated by the scalar field with = 1. This can be an end-point or a future 
attractor (sink) and hence can give the observed late-time accelerated phase. 

P^7 Of these, the points P^ exist only for potentials for which z'^f{z)\z=o = 0. All four 
points are dominated by the kinetic energy of the scalar field with W(p = 1 and = 1 . 
Note that for the points P 2 , the eigenvalue A 3 = Ki given by 

Ki = =fV6{2zf{z) + z^df{z))\z=o- (3.1) 

The nature of the eigenvalues for these points shows that these are the saddle points and 
hence can be associated with a transient phase. 

Ps^s These are matter (barotropic fluid) dominated points. This is reflected in the tCefi 
which is 7 — 1 and = 1- These are also saddle points and hence again a transient 
phase where a phase trajectory can spend some time in the vicinity before moving towards 
a stable attractor. 

P 4^6 These are the radiation dominated points with Wgs = 1/3 and = 1. These are 
also saddle points but given the initial conditions, these can be associated with the past 
attractors (or sources) since our universe was dominated by radiation in the early phase. 
The phase trajectories begin their evolution from this point. 
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Points 

Existence 

Stability 



^0 


Pi 

Always 

z2/(z)|^=o > 0 

-1 

-1 

1 

0 

p± 

?? 

Saddle 

1 

1 

1 

0 

P3 


55 

7-1 

Indet. 

0 

1 

P4 


55 

1 

3 

II 

0 

0 

P5 

55 

55 

7-1 

II 

0 

1 

Pe 

55 

55 

1 

3 

II 

0 

0 

P± 

55 

55 

1 

1 

1 

0 

Ps 

\z*\ < 

See Eq. ( 3 . 2 ) 

S - 3 ) 

1 - 3 ) 

1 

0 

Pg 

> 2 

See Eq. ( 3 . 3 ) 

1 

3 

1 

3 

4 

i? 

0 

Pio 

> \/37 

See Eq. ( 3 . 4 ) 

7-1 

7-1 

37 

1- ^ 

n 


Table 2: Existence and stability conditions for the fixed points for Q = 0. 


Ps-io The properties of these points depend on the potentials since the value of z* is 
determined by the form of f{z). These points can be stable under the following conditions: 


Stability 1 < 7 < | and — -v/Sy < z* < -v/Sy and z*d/* > 0 
condition < 

for Pg I < 7 < 2 and — 2 < z* < 2 and z*d/* > 0 


(3.2) 


Stability condition for Pg s | < 7 < 2 and (z* < —2 or z* > 2) and z*d/* > 0 (3.3) 


Stability condition 4 i < { 9 ^- 2 ) ^ (Stable node) 

for Pro ^ I zj > ^*df* > 0 (Stable spiral). 

(3.4) 

Point Pg is dominated by the scalar field and exists for z^ < 6. Taking 7 = 1 for instance, 
that is, the barotropic fluid is a pressure-less dust, the point is stable for z^ < 3 and will 
give an accelerating phase for z^ < 2 and z*d/* > 0. An interesting aspect of this point is 
that for z^ = 3, the scalar field is able to mimic matter with = WeS = 0. Point Pg is 
usually absent in the analysis done without having a separate radiation fluid component. 
This point offers radiation domination where the scalar field also behaves like radiation 
with W(f) = 1/3. Point Pio is the scaling solution (with Pcp/Pm = const) where neither of 
the components completely dominates the evolutionary dynamics. However, at this point, 
we have = wp , = 7 — 1 . So, an accelerating phase is not possible with 7 > 1. Also, it 
is to be noted that the stability conditions for these points are mutually exclusive so that 
these cannot be stable together unless in a very specific case. 

The above analysis of the critical points shows that we can obtain an accelerated 
expansion provided that the solutions approach the fixed point Pi or Pg in which case the 
final state of the universe is the scalar-field dominated one (H,^ = 1). The scaling solutions 
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Points 

X 

y 

u 

2 : 

Ai 

A 2 

As 

A 4 

Pi 

0 

1 

0 

0 

-2 

-37 

-f±iyi-|zV( 2)|0 

p± 

±1 

0 

0 

0 

3 

1 

3 a T 2/3 

Kl 

P 3 

b 

2/3 

0 

/1 Safe 

V w 

0 

2 

1^2 

- 5 ±^V'c/ 32 -ta 62 

P 4 

20 

3a 

0 

0 

0 

4/32-9a2 

6a 

4/3^—3a& 

6a 

4/32+970 

6a 


P 5 

0 

0 

1 

2a 

0 

2 

-1 

b 

p± 

±1 

0 

0 

2 * 

1 

3 a T 2/3 

6=F\/62;* 

2 


P 7 

/ ^ 

V 3 ^ 

2 

2*\/3 


2 * 

K4 

- 4 z*d/* 

1 1 /I 6 15 

zl 4 

Ps 

2:* 

\/l-# 

0 

2 * 

zS-6 

2 

zi-4 

2 

K 5 

- 2 ^d/* 

Pg 

b 

20 

0 

/ 1 Sab 

V 4/32 

2 * 

80—b\/6zt 

40 

Kq 


Pio 

3a 

0 

0 

2 * 

4/32-9a2 

6a 

4/3^—3afe 

6a 

K7 

Kg 

Pii 

Xu 

2/11 

0 

2 * 

Kg 

Kio 

—szbV s^—At 

4(^6^*-2/3)2 


Table 3: Fixed points and the corresponding eigenvalues for Q = ^j2/3Kl3pm4> where the 
coupling parameter (3 is assumed to be constant. Here o = 2 — 7, 6 = 4 — 87 and c = 47 — 5. 
The point Pn has xn = — 2/3) and yn = V 4/32 + 9a7 — 2\/Q(3z^:I (\/6z* — 2/3). 


Pg and Pio are not viable to explain the late-time acceleration. However, they can be 
used to provide the cosmological evolution in which the scalar field scales proportionally to 
radiation and matter. 

3.2 Conformal coupling between the scalar field and matter fluid 

Having studied the dynamics of the scalar field without any interaction with the matter 
fluid, we introduce a coupling between the two of the form, Q = y^2/3Kl3pm^ in this 
section. This form of the coupling can be motivated [26, 27] by considering a coupled 
gravity and scalar field as a scalar-tensor theory in the Jordan frame and then (conformally) 
transforming to Einstein frame giving a coupled quintessence in normal Einstein gravity. 
In this way the scalar field gets coupled to all matter except radiation for which the stress- 
energy tensor is traceless. The coupling parameter /3 can be a constant or time-dependent 
through an explicit (p dependence. We shall consider a constant /3 for our analysis and 
comment about the variable coupling later. The presence of coupling adds another degree 
of freedom in the parameter space of the quintessence models. Thus, we see the existence 
of new fixed points over and above those present in the analysis without any interaction 
and also modifications to some of the previous points. 

Properties of the fixed points We have points P3, P4, Pg and Pio as the new additions 
and modification of the scaling point Pig (in Table 1) to Point Pn in the new list (Tables 3 
and 4) with the inclusion of coupling. For some of the previous points that still exist with 
coupling, the eigenvalues get modified with the interaction coming in affecting the stability 
conditions with another parameter. 
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Points 

Existence 

Stability 

WeS 




Pi 

V/3 

z'^f{z)\z=o > 0 

-1 

-1 

1 

0 

p± 

V/3 

Saddle/Unstable 

1 

1 

1 

0 

P 3 

Eq. (3.5) 

Always Saddle 

1 

3 

1 

w 

b 

W 

P 4 

|/3| <3(2-7)/2 

Saddle/Unstable 

7-1 + ^ 

1 

4/3^ 

1 _ 4/32 
9a^ 

P 5 

V/3 

Always Saddle 

1 

3 

Indet. 

0 

0 

p± 

V/3 

Saddle/Unstable 

1 

1 

1 

0 

P 7 

> 2 

Eq. (3.10) 

1 

3 

1 

3 

4 

2 ? 

0 

Ps 

\z*\ < a/G 

Eq. (3.11) 

— — 1 

3 

— — 1 

3 

1 

0 

Pg 

Eq. (3.5) 

Eq. (3.6) 

1 

3 

1 

w 

b 

w 

Pio 

|/3| <3(2-7)/2 

Eq. (3.8) 


1 

4/3^ 

9a^ 

~rzrw~ 

9a^ 

Pii 

Eq. (3.13) 

Fig. 1 

Eq. (3.14) 

Eq. (3.15) 



Table 4: Existence and stability conditions for the fixed points with coupling given by 
Q = where the coupling parameter /3 is assumed to be constant. 

Ps^g These are the new points that occur in the presence of coupling between the scalar 
field and matter fluid. The condition for their existence are 

1 < 7 < 4/3 and |/3| > - IO 7 + 8/2. (3.5) 

While P3 is always a saddle point, Pg can be stable under following conditions: 


1 < 7 < 


5 

4 


Stability 
condition < 
for Pg 


f <7< 


4 

3 


' - 10 ^ + 8 

^^ 37 ^ - 10^ + 8 < /j < 

- > - 3V67-476 > Q 

' /? < -^^372 - IO7 + 8 

and and d/* < 0 


(3 > ^^372-107 +8 

and z, > and d/, > 0 

Note the eigenvalues referred to in the Table 3 for these points have 


(3.6) 


^2 = -1/2 


3 (4 - 37 ) 

/3 


{2zf{z) + z^df{z)) 


2=0 


Kq = 


3 (4 - 37)2:gd/^ 
2 f3 


(3.7) 


Both the points have an effective radiation dominated behaviour with tCefr = 1/3 but have 
W(f, = 1. These points are associated with scaling behaviour which can be stable in the case 
of Pg as shown above. 
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P 4 ,io These are also scaling points with = 0 and ~ constant for constant /?. 

The point P 4 is saddle or unstable while Pio can be stable under the following conditions: 




Stability 

and z* < 

condition 1 < 7 < | < 


for Pio 

0</3< 


and z* > 


2\/6/3 


2\/6/3 

The eigenvalues of these points as given in Table 3 are given by 


«^3 = 


2 2P 


3(2-7) 


i2zfiz) + z‘^dfiz)) 


z=0 


97(2 — 7) + 4/3^ — 2 y / dl 3 z ^ 

^7 — -- 77777- 7 -■; ^8 — 


2 2/3z‘^df^ 


(3.8) 


(3.9) 


6 ( 2 - 7 ) ' V 3 2-7 

For certain values of (3, these points can give accelerated expansion with effective equation 
of state parameter given by WeS = 7 — 1 + 4/3^/9(2 — 7 ). 

P 7,8 These points are same as points Ps.g in the analysis without coupling. With the 
inclusion of /3, the stability conditions for these points get modified to 


Stability 
condition < 
for P 7 


4 f /3 < and <z^<-2 and d/* < 0 

^ I /? > and 2 < z^ < and d/* > 0 


^ = 4 
' 3 


I<1<2{ 


P < 0 and z* < —2 and d/* < 0 
/3 > 0 and z* > 2 and d/* > 0 


o ^ _ 97-12 
^ ^ 2^/6 


and z* < —2 or z* > — 


4vT/3 ' 


^ ^ < -2 or z, > 2 


R > 9'7-12 

^ ^ 2v/6 


and 2 :* < — 


97-12 

:* > 2 
or 2 > 2 

9.y_12 U1 2* ^ z, 


> 2;*d/* > 0 

(3.10) 


Stability 
condition < 
for Ps 


/3 < and - 2 < z* < 0 

l<.< 2 andd /.<0 and 


1 < 7 < 2 and d/* > 0 


/3 < ^ 2 ^ 0 < z* < 


13 > and 0 < z* < 2 

The eigenvalues K 4 and for these points are 


/3+V/3W187 
vT 


-4V^/3z* - 97 zJ + 12z2 

K4 = - - n ; K5 = -37 


/3z* + z 2 ; 


(3.11) 


3z2 ’ y 3 

All the properties for these points are same as that of points Pg^g without coupling. 


(3.12) 
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Figure 1: Stability for the point Pn can be assessed graphically from the real part of 
the eigenvalues X± if the other two eigenvalues are negative for the same set of parameter 
values. The plots here show the forking diagrams of the eigenvalues for ranges in < 0 
(left) and 2 ;* > 0 (right) on respectively. The point will be stable for given (3 and 7 for the 
range in z* for which the real parts of X± (the forks here) he below zero. When the real 
part is quasi-constant and un-forked, the eigenvalues are complex. In this case, if the real 
part is negative, the point will be a stable spiral. 


Pii The scaling point Pio gets modified in the presence of coupling and we designate it 
as Pii in the new table. With non-zero /3, the conditions for existence of this point are 


Existence 
condition 
for Pii 


' /3 = 0 

< 0</?< 1(2-7) 
/3>i(2-7) 

\ 


|z*| > a/37 

^ / 3 —// 3^+187 „„ p+yjj 3 ^+l%'y ^ ^ ^ 4/32-972 + I87 

- V6_ a/6 ^ ^ 2V&P 

- a/6 

(3.13) 


for /3 > 0. The signs in the above expressions change accordingly for /3 < 0. The parameters 
rCeff, ^(j, and Vlrn for this point are 


WeE = (7 - 1) + 


2/7 


W<j, 


97 ^ 


1? ^<j) 


v/iz* — 2(3 ’ 

4/3^ -|-187 — 27 / 5 / 32 * 


(3.14) 


-2/32 - 97 -F 7/6/32;* ' (^2* - 2 / 3 ) 

The complicated nature of the eigenvalues: 

3/3 ,3/ ^ 37/67 2 ^d/ 


2 -= (3.15) 


Kg = 7 


7 / 62 * - 2/3 2 


+ o “ 2; Kio — — 


7 / 62 ;* — 2/3 


; A± = 


—s ± — 4t 

4 ( 7 / 6 ^* - 2/3)2 


(3.16) 
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where s and t are given by 


s = 18az^ + 24/3^ + — 24\/6/3z* 


(3.17) 


and 

t = - 1296a^/3%2 - 'dl2a^-i^zl + 1296\/6a2/372z* - 576a/3S - + 12960/3^7^ 

+ 548072 ;^ — 432\/6o/37z^ — 864a/3^7z^ — 97207^2^ + 72O\/6a;0^7z* + 1296\/6a/37^2;* 
- 1296\/6o/3722* - 576/3^7^ + 576/3^7 - \^^^/l^zl + 1152/3^4 - bl&y/Q^^zl 
+ 21QVQzl + 7Q8l3^zl - 2160/3%2z2 + 172815‘^'yzl - 64\/6/3®z* + 864\/6/3^72z* 
- 864 V 6 / 337 Z* (3.18) 

makes it difficult to assess the conditions for stability for this point analytically. The point 
will be stable node if all the eigenvalues are negative and a stable spiral if \± are complex 
with the real part being negative. We depict this for some of the cases in the Fig. 1. 
The plots show the forking diagrams of the eigenvalues for ranges in 2 ;* < 0 (left) and 
z* > 0 (right) on respectively. The point will be stable for given /3 and 7 for the range 
in z* for which the real parts of \± (the forks here) lie below zero. When the real part is 
quasi-constant and un-forked, the eigenvalues are complex. In this case, if the real part is 
negative, the point will be a stable spiral. 

We assumed the coupling parameter /3 to be constant but generally it is (j) dependent 
or time-dependent. In this case, we need another equation to close the autonomous system. 
Defining another variable s = Kcj), this gives ds/dN = \^x. For the fixed points of the new 
system, we then need the right hand side of this equation to vanish as well. This makes 
a: = 0 necessary for all the fixed points. Then, the subset of fixed points listed in Table 3 
for which x = 0 form the instantaneous fixed points for the system with changing /3. The 
exact analysis would then require the functional form of the variable coupling to assess the 
exact dynamics. 

Interpretation of the fixed-point analysis. Before we proceed further, we shall com¬ 
ment briefly about the methodology used and the results obtained in this section. As noted 
earlier in the previous sections, the cosmological system presents a complicated non-linear 
problem in which the dynamics of the scale factor depends on the constituent fluids in 
a non-trivial manner. The inclusion of a scalar field to provide a dynamical dark energy 
constituent further adds on to the complexity of the problem. Then, to ascertain the dy¬ 
namics of the system, we need to solve simultaneously, the Friedmann equation and the 
Klein-Gordon equation of the scalar field which gives the evolution of the scale factor and 
the scalar field given the initial conditions and the potential of the scalar field. The system 
is thus far from being solvable analytically, in general, and we need to resort to numerical 
techniques. 

However, the otherwise analytically opaque dynamics can still be understood, albeit, 
in the manner of computing the fixed points of the non-linear equations. This is a fairly 
standard technique in non-linear dynamics which has been employed in dark energy models 
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as well (see, for example, ref. [9]). The fixed points are to be taken as the stages a non¬ 
linear system may go through during its evolution given the initial conditions. We adopted 
a similar procedure in this section with some differences and additions, viz., inclusion of 
radiation fluid, keeping the matter fluid as a barotropic one (with non-zero pressure), form 
of the potential for a scalar field to be generic and, including a coupling between the 
scalar field and matter sector. This, in a way, stretches the scheme of things to be as 
generic and tractable as possible. Thus, the fixed points obtained in this section form a 
super-set of all possible stages that a (/>CDM model (with and without coupling) can go 
through during its evolution. Then, for a specific potential of the scalar field, the allowed 
set of fixed points for that potential is a subset of this super-set and these fixed points 
give an instantaneous behaviour of the system in the phase space. Another aspect to be 
understood is that, although given a potential we can have a set of fixed points for the 
same, the points a system will traverse on its trajectory in the phase space depends on the 
choice of the initial conditions. Thus, the set of fixed points and their nature gives us a 
possible set of trajectories in the phase space which in itself is worth knowing, given that we 
can hardly infer anything analytically otherwise. The estimation of parameters and testing 
the concordance of models with the observations, in the end, is a numerical exercise. 

4 A numerical peek at the cosmological dynamics 

Having analysed the radiation-matter-quintessence dynamical system using the linear sta¬ 
bility analysis, we now look at the exact (using numerics) dynamics for a few examples. We 
take the function f{z) to be of the form: z and — z)lz^ with and without cou¬ 

pling and solve the dynamical equations for the observed cosmological parameters. These 
pertain to the potentials H(<?i) of the form: Vq exp(—/uj!)^/ 2 Mp), Vb exp(exp(—/r(/)/Mp)) and 
Vq exp(— exp{(j)/Mp) — ^(p/Mp) respectively. The first two cases do not have a finite z* and 
the third case falls under the category z* = /r with z'^f{z)\z=o 7 ^ 0. Finally, we also consider 
a case f{z) = {^ — z) which has z^ = fi and z'^f(z)\z=o = 0. This suffices to illustrate and 
explain the behaviour of models in accordance with the fixed point analysis of the previous 
section. We also assume that the barotropic fluid is perfect matter, that is, we take 7 = 1 
for all the analysis henceforth. 

We show the evolution of density parameters (H^, and density of the com¬ 
ponents, equation of state parameters and rCeff for the three cases [i\jz^ ^ 1 ^ 2 1z and 
(^3 — z)lz^ in Figure 2 with and without coupling. We also give the evolution of the co¬ 
moving Hubble radius, the evolution of the scalar field and trajectories of the system in the 
phase space in Figure 3. The parameters are taken to be /ii = 1, /r 2 = 39 and /is = 100 
and the coupling constant j3 = 0.1 assuming a weak coupling. We evolve the system by 
setting the initial conditions at the epoch of nucleosynthesis which corresponds to A ~ —22 
(although the plots for the Hs are clipped to show from A ~ —10 since initially it is all radi¬ 
ation!) from temperature considerations. For the initial conditions, we choose = 0.9 and 
7 (^. =4/3 and an initial slope parameter Zi so that the scalar field mimics radiation at the 
nucleosynthesis epoch. For determining initial normalized Hubble parameter hi = Hi/H q 
and initial we use computed from Tcmb = 2.725 K and take = 0.286 [44]. 
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N 


Figure 2: Numerical peek at the cosmological dynamics I. The plots show the 
evolution of the density parameters (top row and middle-left) for f{z) given by z, ^ 2 ! 
and ( 1 U 3 — z)lz^ with (dashed) and without (regular) coupling. The parameters are /ii = 1, 
/X 2 = 39 and = 100 and the coupling constant /3 = 0.1. All three cases show tracker 
behaviour for the scalar field. This is particularly evident in the plot showing the evolution 
of density (middle-right) for f{z) = (ps — z)/z^. Finally, we show the evolution of and 
WeS (bottom row) for the three cases. 


With this choice of initial conditions, we see that for all three cases, f{z) = fii/z'^, ^^ 2 !z and 
(/ra — z)lz^, the dynamics of the scalar field is of a tracker type. This is particularly evident 
in the evolution of density in Figure 2 (middle-right) for f{z) = (//s — z)lz^. The density of 
the scalar field falls closely with radiation till it dominates and then scales with matter after 
radiation-matter equality and finally exits to give the late-time acceleration. The comoving 
Hubble radius (top rows in Figure 3) starts off as a line of unit slope featuring the radiation 
era which goes into the matter phase and then tilts down again with a unit slope implying 
accelerated expansion. We also show the evolution of w^j, and foi' the three cases. The 
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Figure 3: Numerical peek at the cosmological dynamics II. The plots show the 
evolution of the comoving Hubble radius with (top-right) and without (top-left) coupling, 
the evolution of the scalar field (bottom-left) and trajectories in the phase space (bottom- 
right) with (dashed) and without (regular) coupling. 


effect of coupling is seen to be more pronounced on than any other quantity. Such an 
effect was also shown in ref. [45] in the case of specific potentials. It is interesting to note 
therefore if we can have any observables which are local in redshift, it can serve to constrain 
the model, f(z) and the coupling. Unfortunately, the current observables like luminosity 
distances integrate out these features leading to degeneracy between many models. 

Finally, we consider an example with f(z) = ^ — z without coupling. For this case, 
noted earlier as well, we have z^: = and z'^f{z)\z=o = 0. We take different values for the 
parameter = 2, \/3, y/2, 1 and —1 to illustrate different aspects of evolution of the system 
as per the stability of the fixed points. We plot the phase space variables (x, y, z, u) and 
rCeft for these cases in Figure 4. The plots confirm the following aspects: 


1. For y = 2, it is seen from the conditions of stability, that only point Pio in Table 1 
is stable. Hence, we see in Figure 4 (Top-left), the phase space variables (regular 
lines) approach the values for Pio with x = y ~ 0.61, u ~ 0 and z = 2 and tUgg 
(Bottom-right) goes to zero at late-times. 
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Figure 4: The plots in this hgure display the evolution of phase space variables and tCefr 
for f{z) = — z without coupling. We take different values for the parameter fj, = 2, \/3, 

\/2, 1 and —1 to illustrate different aspects of evolution of the system as per the stability 
of the fixed points. The plots show and confirm the scaling behaviour, non-accelerated and 
accelerated phases and the de Sitter phase in the late-time evolution for these values of 
parameter /x which is the zero of f{z). 


2. For /X = y/3, the system evolves to point Pg in Table 1 which is stable in this case and 
the scalar held mimics matter. The phase space variables (dashed lines) in Figure 4 
(Top-left) approach the values for Pg with x = y k, 0.707, u « 0, z = v/S and WeE ~ 0 
(Bottom-right) at late-times. There is no late-time acceleration in this case. 

3. For y = \/2, the stable point Pg shows critical behaviour as WeS at late- 

times. The phase space variables (regular lines in Figure 4, Top-right) approach the 
values X = 0.577, y = 0.816, tt ps 0 and z = y/2. 

4. Point Pg shows accelerated behaviour for /x = 1. The phase space variables (dashed 
lines in Figure 4, Top-right) approach the values x = 0.408, y = 0.912, u ~ 0 and 
z = 1 while xCeff —>■ —2/3 at late-times. 

5. Flipping the sign of y with /x = — 1, point Pi in Table 1 becomes the stable node. 
Hence, we see in Figure 4 (Bottom-left), the phase space variables approach the values 
x = 0 , y = l, u = 0 and z = 0 asymptotically with tCefi — )• — 1 giving asymptotic 
de Sitter behaviour. This can be tuned to give the correct cosmological parameters 
at the present epoch but in the ‘thawing’ class of quintessence models with a shallow 
slope. 
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The numerical computation of the dynamical equations thus validates the fixed point 
analysis to ascertain the background evolution of the quintessence models. The potential 
for the scalar field affects the background dynamics through the function, f{z). It would 
be interesting to see if there is a way to reconstruct this function given the observational 
data. 


5 Summary and Discussion 

‘Quintessence’ models explore the idea of using a scalar field as the dark energy component 
in the cosmological soup. The extra scalar degree of freedom also arises naturally if one 
considers modifications to general relativity in various ways and the resulting models can 
also pick up a coupling between the scalar field and matter fluid leading to a coupled dark 
sector dynamics. Between all this, the potential for the scalar field is motivated either 
phenomenologically or from higher dimensional theories. The question we ask ourselves is, 
“Can we understand the background cosmological dynamics with a generic potential for the 
scalar field?”. 

It turns out that it is indeed possible to do so in a way that can encompass a large 
class of potentials. The cosmological dynamical equations can be transformed into an 
equivalent autonomous system of differential equations. The non-linear system can be 
analysed analytically only through the linear stability analysis which includes computing 
the fixed points for the system and assess their stability. Even at this level, it is seen to be 
quite an effective tool for understanding background dynamics since these fixed points can 
be associated with various cosmological solutions that can decide the fate of the evolutionary 
dynamics. 

We performed such an analysis by considering a three eomponent system with a scalar 
field, matter and radiation fluids but keeping the potential for the scalar field generic. 
We also considered a coupling between the scalar field and matter fluid of the form Q = 
which can be motivated by considering a non-minimal coupling in the Jordan 
frame and then transforming the action to Einstein frame. It turns out that the equivalent 
dynamical system depends only on a function f(z) of the relative slope of the potential 
which can be classified into three main categories for a large class of potentials beyond 
simple exponential type. 

We find that the fixed points for the three fluid system describe different cosmological 
solutions and can be associated with source/sink/transient phases of evolution depending 
on their conditions of stability. We see the existence of de Sitter dominated point, scaling 
solutions in which the scalar field behaves like radiation or matter thus leading to tracker 
type behaviour, matter/radiation dominated points. These points have different stability 
conditions and cosmological parameters. For some of the cases, the de Sitter point can 
become a global attractor and can be associated with the late-time acceleration. The effect 
of coupling renders few more scaling points and also modifications to the stability conditions 
of the persisting points. The coupling is introduced to account for the coincidence problem 
of ACDM model and allows for the interaction and transfer of energy between the dark 
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matter and quintessence field. The analysis, to be stressed again, is done independent of 
the form of the potential for the scalar field. 

Finally, we used numerical techniques to illustrate and confirm the analysis by taking 
a few examples. We considered examples from each of the cases from the classification by 
taking the function f{z) of the form (/^ ~ (/^ ~ ^') with and without 

coupling. For the choice of hand-picked initial conditions, we see the realistic evolution 
of the cosmological parameters for these cases. The first three examples show tracker 
type behaviour for the scalar field while the last one can be realised with the observed 
cosmological evolution only under thawing class of scalar field models. Further, within 
the last case, we illustrate the aspects of the fixed point analysis in different settings. We 
observe that the coupling shows major effect on the equation of state parameter of the 
scalar field and also allows for smoothening of the density evolution for the scalar field 
in the scaling regime. However, most of the current background-level cosmological probes 
integrate out these effects along the redshift leading to degeneracy between many models. 
It would be interesting to see how we can employ the parameter-estimation methods and 
study the perturbed evolution also in a generic way. This will have two-fold motive: (i) to 
study deviations from ACDM and (ii) if possible, reconstruct the potential or the function 
f{z) for the quintessence field. 
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